Dielectronic Resonance Method for Measuring Isotope Shifts 
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Longstanding problems in the comparison of very accurate hyperfine-shift measurements to the- 
ory were partly overcome by precise measurements on few-electron highly-charged ions. Still the 
agreement between theory and experiment is unsatisfactory. In this paper, we present a radically 
new way of precisely measuring hyperfine shifts, and demonstrate its effectiveness in the case of the 



hyperfine shift of 4si/2 ^nd ^Pi/2 



/p^i, 



It is based on the precise detection of dielectronic 



resonances that occur in electron-ion recombination at very low energy. This allows us to determine 
the hyperfine constant to around 0.6 meV accuracy which is on the order of 10%. 

PACS numbers: 31.30.Gs, 34.80.Lx 
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The hyperfine structure (hfs) in atoms reflects fun- 
damental high-order effects in the interaction between 
the bound electrons and the nucleus. The energy shift 
contains contributions from Quantum Electro-Dynamics 
(QED) and from the nuclear charge and magnetic mo- 
ment distribution 1]. These effects have been studied 
over decades by spectroscopy, mainly on neutral atoms 
or singly charged ions. In such many-electron systems, 
the atomic structure is complex. The theoretical pre- 
dictions are further hampered by imprecise knowledge 
of nuclear magnetic moments which, if measured in a 
magnetic field, is shielded due to diamagnetic effects [l|, 
which can only be estimated and corrected for through 
calculations. Only recently, experimental approaches to 
measure the hyperfine structure of heavy highly-charged 
systems were opened y, y, Ij, |5|, giving a new way 
to obtain nuclear moments and to test QED in a new 
regime. Measurements were done for several hydrogen- 
like ions: 2093^82+^ ^e5^o''e+ ^ I83,i87j^e74+^ 203,205^180+, 

and 2°^Pb®^+|l|l|ll3. In all cases the theory and ex- 
periment do not agree in a satisfactory way. The largest 
theoretical uncertainty comes from the unknown distri- 
bution of nuclear magnetic moment (Bohr-Weisskopf ef- 
fect), which is complicated to predict due to nuclear po- 
larization [g, but also precise knowledge of the charge 
radius is crucial. On the other hand, a determination 
of nuclear magnetic-moment distribution is impossible 
due to uncertainties in the QED prediction. Although 
the measurements with H-like ions are quite accurate 
(2 X 10"'' for 2"7pb 0), it is still essential to find al- 
ternative approaches to measure hyperfine structure for 
different ion charges where the above mentioned effects 
contribute differently. 

Here we demonstrate a first measurement of the hy- 
perfine structure in ^o^Pb by observing dielectronic re- 
combination (DR) resonances for both 208p|-,53-i- j^^ero 
nuclear spin) and ^o^pi^ + (nuclear spin one-half). Us- 
ing the two Pb isotopes in charge state 53+, one clearly 
observes hyperfine effects in 4si/2 and 4pi/2 F=l and 



F=0 states. The importance of this observation is that it 
provides means to easily measure the hyperfine interac- 
tion of 2*''''Pb ~'', and other ions with a different number 
of electrons, as compared to what is currently possible. 
By measuring the same nucleus with different numbers 
of electrons one can measure their influence on the mag- 
netic interaction and check the theories including QED 
and probably do even better by extrapolating to low elec- 
tron numbers. A similar idea, combining the Is hyperfine 
shift in H-like and the 2s shift in Li-likc ions to extract 
QED contributions, was recently proposed in a theoret- 
ical paper lZ|. A further application of this method can 
be found in the new accumulator rings for radioactive 
nuclei to measure their nuclear magnetic properties by 
low-energy electron scattering in the cooler 8J| . 

The spectroscopic principle of using DR resonances is 
that the valence electron is excited by the capture of a 
free electron into a Rydberg state. The measured reso- 
nances are near threshold, thus the Rydberg energy must 
balance to a large part the excitation energy of the va- 
lence electron. Higher order effects in the interaction of 
the Rydberg electron with the core electrons are weak 
and can be estimated to excellent accuracy. The mea- 
surement is checked by comparing the DR resonances 
with those of ^o^Pb ^ where no hyperfine splitting ex- 
ists. In the 208p|-, + ^g^gg ^Q reached before an accuracy 
that enabled us to obtain spectroscopic data for stringent 
tests of many-electron QED effects in advanced atomic 
structure calculations 10]. 

Selected ^'''^Pb or ^o^pb isotopes are injected as singly 
charged ions from an isotope separator, at the Manne 
Siegbahn Laboratory in Stockholm, into an electron- 
beam ion source, where the highly charged Pb ions are 
created. From there they are injected into the heavy- 
ion storage ring CRYRING, and accelerated to several 
MeV/amu. After electron cooling the ion beam for 1 
s, the recombination rates as a function of relative en- 
ergy between electrons and ions are measured. This is 
achieved by ramping the electron energy up and down. 



from the cooling condition, to a maximum energy in the 
center-of-mass frame so that recombination spectra are 
obtained both with the electrons being slower and faster 
than the ions. Thus, we can check the space charge of 
the electron beam and the longitudinal drag force ex- 
erted on the ions by the electrons during detuning and 
make appropriate corrections |llLll2| . These corrections 
are important in order to obtain accurate energy scales 
in the spectra. The energy calibration was checked thor- 
oughly and found to be accurate to within 0.1 meV at 
energies < 0.3 eV. For these low energies one takes full 
advantage of the small energy spread of the adiabatically 
expanded electron beam jl3| and obtains a resolution on 
the order of 1 meV. From fits to the resonances, we con- 
firm the longitudinal and transversal temperatures of the 
electron beam to be kT; =0.1 mcV and kTf — 2 meV, 
respectively. 

The consistency of the measurement is checked by 
switching between the 207p|353-H ^^^ 208p|353-(- 1,^^^^^^ 

Small changes in the magnet field strengths were made 
for the two isotope ions to have the same velocity (deter- 
mined by electron cooling) and the same circulation fre- 
quency in the ring (measured by the Schottky frequency) . 
We also compare the 208pj-,53-i- ^^^^g^ measured in this run 
with those measured earlier (with kT; = 0.08 meV and 
kTt = 1 meV) and published in Ref. [I^ (see Fig. ^. 
Considering the small temperature difference, these data 
sets agree within the statistical error (its size is seen in the 
spread of the points). In Fig. ^we also show the recom- 
bination rate coefficient that was obtained with ^°^Pb ^ 
stored in the ring. The difference between the two sets 
of experimental spectra must then be due to the isotope 
effect, in particular the hyperfine shifts of the F — and 
1 states in Asi/2 and 4pi/2 of ^°'^Ph^^^ . 

We have evaluated the isotope dependent contribu- 
tions to Pb^^+ (4^1/2) using both Relativistic Many- 
Body Perturbation Theory (RMBPT) ^ and a Multi- 
Configuration Dirac-Fock (MCDF) computer pack- 
age [ij, [13 . We have carefully checked the agreement 
between the methods. The dominating effect, listed in 
detail on the first six lines of Table Ql is due to the cou- 
pling between the total electronic angular momentum 
and the nuclear magnetic moment of the spin one half 
^°^Pb nucleus. Most of the contribution is accounted for 
already with a single-configuration DF calculation using 
a point nucleus, as seen on the first line. The lead wave 
function has then been evaluated in a potential from an 
extended nucleus, which results in a ^ 10 % change in 
the 4s hyperfine structure, as seen on the second line 
of Table Q] We use a Fermi distribution of the nuclear 
charge with thickness parameter 2.3 fm and a half den- 
sity radius adjusted to produce a charge root mean square 
radius (rms) of 5.4942 (0.0013) fm [l3|. The wave func- 
tion has also been allowed to adjust to the Breit interac- 
tion between the electrons and to the Uehling potential 
(responsible for the bulk of the vacuum polarization). 



but the effect on the hyperfine structure is very small. 
All these contributions agree within the displayed digits 
between the two computational methods. The polariza- 
tion of the closed shell core by the valence electron give 
rise to important corrections to the hyperfine structure. 
These effects can only be accounted for when spherical 
symmetry is abandoned, but they are still not true corre- 
lation effects. We have calculated this core-polarization 
using RMBPT(fourth line). The dominating part of the 
core-polarization is accounted for already in second or- 
der (one order in the hyperfine structure and one order in 
the Coulomb interaction). The result change with only 
a few /ieV when higher order core polarization is taken 
into account. Since also true correlation, calculated with 
either RMBPT or MCDF, gives negligible contributions, 
see Table U it is evident that the many-body aspect of 
p|-,53-(- jQgg j^Q^ limit the potential accuracy of the de- 
termination of the hyperfine structure. An approximate 
Bohr-Weisskopf effect is included, assuming the same 
mean spherical radius for the magnetic moment distri- 
bution as was used for the charge distribution, following 
e.g. Ref. 'ISl. Also here the two computational methods 
agree to within one unit in the last quoted digit. 

An additional isotope dependent contribution origi- 
nate from the slightly different potentials produced by 
nuclei with different rms (volume shift). According to 
Angeh [13 the rms of ^"^pb is 0.0068 (0.0002) fm smaller 
than that of ^°*Pb. This leads to ^ 1.6 meV larger bind- 
ing energy for the 4s state of ^"^Pb, corresponding to 
- 0.2 meV /lO"^ fm, as seen in TableU 

The experimental evaluations of the magnetic moment 
of e.g., ^°^Pb have been critically evaluated in Ref.Jj, in- 
dicating that the true value should be within error bars of 
the Atomic Beam Magnetic Resonance (ABMR) result by 
Brenner ^/ = 0.5918 (14)^Ar ilg| and the Nuclear Mag- 
netic Resonance (NMR) result by Lutz and Strieker |19l |. 
quoted in Ref. 1] as /i/ = 0.5925 (6)/i,w. While the optical 
pumping result by Gibbs and White [2|, which is « 1.3% 
below Ref. jl8| . seems to suffer from s yste matic errors. 
Here we have used the value from Ref. [isf. 

As in ^"^Pb |1D|| the recombination goes through 
resonances in Pb 
Pb52+ 



52^ 



Pb53+ (3di04si/2) + e ^ 
(3d^°4pi/2l8£j) ,. The term splitting of these ex- 
cited states range from a few meV for the states around 
the ionization threshold, j = 21/2, of 208pb^^+[l3 to 
practically zero for the highest angular momenta. When 
considering hyperfine structure we note first that the J- 
dependent term splitting and the hyperfine splitting of 
the 4j5i/2 orbital energy are of the same order of magni- 
tude. The energy levels are then expected to be governed 
by intermediate coupling and will not scale linearly with 
the hyperfine shift; each level will be an individual mix 
of the two hyperfine components of the inner electron or- 
bital. The hyperfine interaction involving the outer 18£j 
electron will be much smaller due to the 1/n^ scaling of 
the radial matrix elements and will be neglected in the 



following. Starting from the f4pi /^IS ^i) , energy levels 
in the absence of hyperfine interaction [1^ , using the as- 
sumption that only terms (J) belonging to the same con- 
figuration will mix due to the hyperfine interaction, and 
using the knowledge of the splitting of the 4pi/2 iu the 
absence of the Rydberg electron, it is straight forward to 
generate the hyperfine perturbed resonance energies with 
standard angular momentum theory. Without hyper- 
fine interaction each 4pi/2l8£j configuration forms two 
J-terms. With hyperfine interaction we get four energy 
levels of which two mix the J-terms. The recombination 
cross section can now be constructed from the energy 
positions relative the initial state of Pb^^^ , the autoion- 
ization rates and the stabilizing radiative rates to bound 
states of Pb^^"''. The autoionization rate varies strongly 
and is governed by the actual mixture just discussed. The 
radiative rate is for these high angular momenta states 
dominated by the probability for the 4p — > 4s transition 
and is thus slowly varying and nearly constant for all the 
investigated resonances. The position of a resonances is 
affected by the positions of the hyperfine perturbed ex- 
cited states as well as of the shifted initial state. The 
ground state of Pb^^^4s is with hyperfine interaction 
split into one state with F=0 and one with F=l. The 
F=l state has higher energy, but its decay to F=0 is 
very slow. The expected life time is 2.3 x 10"* s y] which 
is much longer than the time from ions production to the 
measurement (~ 5s). We have thus assumed a statisti- 
cal population of the two Pb^^^ (4si/2)^ states and the 
recombination cross section from the two possible initial 
states have been calculated and added with statistical 
weights. This cross section has subsequently been folded 
with the electron beam velocity distribution mentioned 
above. 

The nuclear spin dependent part of the hyperfine en- 
ergy shift is A^^^S = a/2[F(F-M)-/(/+l)-J(J+l)], 
where a is the hyperfine splitting constant for a given 
state, J is the electronic angular momentum, / is the nu- 
clear spin, and F is the total angular momentum. Using 
the fact that Sa^p « 045 to within 2 — 3%, as shown in Ta- 
ble ^ we generate a synthetic spectrum derived from the 
resonance parameters as described above, and test the 
sensitivity with respect to small variations in a = a^p. 
The theoretical shifts in Table HI are approximately re- 
produced by a = 4.8 meV. Fig 12 shows the spectrum 
generated with the theoretical shifts as well as with a val- 
ues deviating from it with around ±0.7 meV. The spec- 
trum generated from the calculated shifts seems to fit the 
data best. Note that different peaks in the spectrum are 
rather differently affected when a is changed, refiecting 
the non- linear shift of the resonances. With a = 5.5 meV 
the synthetic spectrum shows much more pronounced res- 
onances and with a — 4.0 meV the dip around 0.01 eV 
has nearly disappeared. We estimate from this procedure 
that with the present energy spread of the electron beam 
we can get an accuracy in determining a, from a fit to the 



resonances, to the order of ±0.3 meV. Here the statistical 
error was taken into account. With a specially designed 
electron target of lower energy spread, one could reduce 
the uncertainty further. Fig 13 shows the positions of all 
the resonances and the cross section before folding with 
the electron temperature. 

At the present stage we are also limited in our accuracy 
by the uncertainty in describing the resonances. This 
problem already exists in the ■^°^Pb '^ case [l(]j. There 
is a so far unexplained discrepancy between theory and 
measured rate coefficient at energies below 1 meV. That 
may cause the systematic uncertainty to reach 1 meV. 
Several resonances above 1 meV are well described by 
the calculated series. Another, unexplained discrepancy 
appears at energies larger then 0.03 eV (see Fig. 3). As 
can be seen from Fig. 3 the highest energy resonances 
contribute insignificantly to the cross section. The rea- 
son might be that the autoionization rates approach zero 
for these high angular momenta states, but even small 
electromagnetic fields in the cooler could change this by 
inducing a mixture with states of low angular momen- 
tum. This possibility has to be investigated further, but 
will not affect the conclusions here since in this range one 
is not very sensitive to a hyperfine splitting in the order 
of mcV. 

In conclusion, the hyperfine splittings in 4pi/2 and 
4si/2 state of the Cu-like ^"^Pb isotope were measured by 
a new method to an accuracy of the order 10% of the hfs 
constant. Instead of the generally applied optical spec- 
troscopy, the splitting was observed in dielectronic res- 
onances in low-energy electron-ion recombination. This 
opens new possibilities to measure this important physi- 
cal quantity in a new energy range and ionization stage 
for determining the different contribution to its value. 
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FIG. 1: Measured dielectronic recombination spectra with 
20^Pb5^+ and with 208 Pb53+. 
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FIG. 2: Calculated recombination cross sections for Pb 
with different values of the hyperfine constant, a, as explained 
in the text. 
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FIG. 3: Calculated resonance positions for initial states 
^"'^Pb^^^ (4S1/2) „_„, upper sequence of levels, and 



207 



F=0' 



Pb^^"*" (4S1/2) „_, , lower sequence of levels. The grey 
line shows the cross section assuming statistical population 
of these two hyperfine structure components. From this cross 
section the rate coefficient, the black solid line, is generated. 
It should be compared to the measured rate coefficient for 
207p|^53+ gjjo^jj as black dots. 



TABLE I: Binding energy of ^o'^ph^'i+ (4Ji/2)^ relative the 
(4ii/2) states in 208PI353+ (jjjeV). 

4Si/2 4pi/2 

F^O F^l F=0 F=l 

Contributions ~ ^j" 

point size fij 

Dirac-Fock point nucleus -12.16 4.05 -3.66 1.22 

A nuclear size 1.37 -0.46 0.12 -0.04 

A Breit and Vac. pol. -0.03 0.01 0.03 -0.01 

A core-polarization' -0.34 0.11 -0.12 0.04 

A correlation -0.03 0.01 0.04 -0.01 

distributed fij 

(Bohr- Weisskopf effect) 0.23 -0.08 0.02 -0.01 

Sum -10.97 3.66 -3.57 1.19 

Contributions independent of /i/ 

Volume shift -1.59 -1.59 -0.08 -0.08 

Total -12.56 2.06 -3.65 1.11 

''iJ.i/lJ.N = 0.5918, Ref. [13 

''Includes polarization due to both the Coulomb and the Breit 
interaction, the latter contributes with 10-14% 



